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Abstract. Recent results on the fluid dynamic limits of the Boltzmann equation based on the DiPerna-Lions theory of 
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In memory of Carlo Cercignani (1939-2010) 

Relating the kinetic theory of gases to their description by the equations of continuum mechanics is a question that 
finds its origins in the work of Maxwell [30]. It was subsequently formulated by Hilbert as a mathematical problem 
— specifically, as an example of his 6th problem on the axiomatization of physics [21]. In Hilbert's own words 
"Boltzmann's work on the principles of mechanics suggests the problem of developing mathematically the limiting 
processes which lead from the atomistic view to the laws of motion of continua". Hilbert himself studied this problem 
in [22] as an application of his theory of integral equations. The present paper reviews recent progress on this problem 
in the past 10 years as a consequence of the DiPerna-Lions global existence and stability theory [12] for solutions 
of the Boltzmann equation. This Harold Grad lecture is dedicated to the memory of Carlo Cercignani, who gave the 
first Harold Grad lecture in the 17th Rarefied Gas Dynamics Symposium, in Aachen (1990), in recognition of his 
outstanding influence on the mathematical analysis of the Boltzmann equation in the past 40 years. 

THE BOLTZMANN EQUATION: FORMAL STRUCTURE 

In kinetic theory, the state of a monatomic gas at time t and position x is described by its velocity distribution function 
F = F(t,x,v) > 0. It satisfies the Boltzmann equation 

d,F + vV x F = tf(F), 

where c tf(F)(t ,x, v) := ^(F(t,x, -))(v) is the Boltzmann collision integral defined for each continuous, rapidly decay- 
ing function / = /(v) by 

:= / / ,(/(vO/K)-/(v)/(v,))£ |(v-v*)-<o|dv*d<0, 

JjR 3 xS 2 

assuming that gas molecules behave as perfectly elastic hard spheres of diameter d. In this formula, we have denoted 
v = v (v, v*,fi>) := v — (v — v*) • coco, = v„(v,v*,Ci)) := v*+ (v — v*) • coco, |©| = 1. (1) 

,2 

Molecular interactions more general than hard sphere collisions can be considered by replacing %-|(v — v*) • CO\ with 
appropriate collision kernels of the form b(\v — v*|, | ^Zl*\ ' <*>!)■ I n m i s P a P er ' we restrict our attention to the case of 
hard sphere collisions to avoid dealing with more technical conditions on the collision kernel. 

Properties of the collision integral 

While the collision integral is a fairly intricate mathematical expression, the formulas (1) entail remarkable symme- 
try properties. As a result, the collision integral satisfies, for each continuous, rapidly decaying / = /(v), the identities 

/ ^(/)dv = 0, / V(f)v k dv = 0, k= 1,2,3, and / tf(f)\v\ 2 dv = 0. (2) 

JR 3 JR 3 JR 3 



The first relation expresses the conservation of mass (or equivalently, of the number of particles) by the collision 
process, while the second and the third express the conservation of momentum and energy respectively. 

Perhaps the most important result on the collision integral is Boltzmann 's H Theorem: for each continuous, rapidly 
decaying / = /(v) > such that In / has polynomial growth as |v| — > +°°, 

/ %?(/)ln/dv<0, and / %?(/)ln/dv = <=> <*f(/) = <=> / is a Maxwellian, (3) 

JR 3 JR 3 

i.e. there exists p, 9 > and u e R 3 such that 



Dimensionless variables 

Fluid dynamic limits are obtained as properties of solutions of the Boltzmann equation under appropriate scaling as- 
sumptions. We therefore recast the Boltzmann equation in dimensionless variables, so as to identify the dimensionless 
parameters that control the scalings of the time and space variables, following [4, 36]. 

First we choose a macroscopic length scale L (for instance the size of the container where the gas is enclosed, or of 
an object immersed in the fluid, or the typical length scale on which the variation of macroscopic fluid quantities is 
observed), as well as a macroscopic observation time scale T (i.e. the time scale on which the evolution of the fluid 
quantities is observed.) 

We next define reference scales of density p and temperature 9 by setting 

jjFdxdv = pL\ JJvF6xdv = 0, jj \\v\ 2 FAxAv = \ p6 . 
The collision time scale T c is defined in turn by 

T /// ^(p A e)M-^(p,o,e)( v *)K v_v *)' fi) l dvdv * dfi);= y> 

while the acoustic time scale is defined by T a = L/Vd. The dimensionless variables are i = t/T , x = x/L, and 

v = v/V?, while the dimensionless distribution function is F = 9 i ^ 2 F /p. 

Introducing two dimensionless parameters, the Strouhal number Sh = T a /T and the Knudsen number Kn = 

CT c /T a =2/d 2 Lp with 

c^^yyiv-v,| e -(H 2 +iv.i 2 )/2 dvdv ^ 

we see that the Boltzmann equation in dimensionless variables takes the form 

Sha f F + v-V je F = ^(F), (5) 
where the dimensionless collision integral is 

e £{P){i,x,v) = [[ (F(i,x,v)F(i,x,v[)-F(i,x,v)F(i,x,v^)\(v (6) 

JjR 3 xS 2 

Obviously, the dimensionless collision integral ^(F) satisfies exactly the same properties as the original expression 
ff(F), i.e. the conservation laws of mass, momentum and energy (2) and Boltzmann's H Theorem (3). 

Henceforth, we always consider the Boltzmann equation (5) in dimensionless variables, dropping all hats for 
notational simplicity. Thus, the conservation properties of the collision operator (2) imply that rapidly decaying (in v) 



solutions of the dimensionless Boltzmann equation (5) satisfy the following local conservation laws: 

Shd, / Fdv + div x / vFdv = (conservation of mass), 

Sh<9 ( / vFdv + div x / v<givFdv — (conservation of momentum), 

Sh(9 ( / l\v\ 2 F dv + div x vj\v\ 2 Fdv~0 (conservation of energy). 
JR 3 Jr3 

Likewise, Boltzmann's H Theorem implies that solutions F > of the Boltzmann equation that are rapidly decaying 
while InF has polynomial growth as |v| — > +°° satisfy the differential inequality 



(7) 



Shd, / F\nFdv + div x [ vFlnFdv=— [ %?(/)ln/dv < 0. 
JR3 JR3 Kn Jr3 



(8) 



THE COMPRESSIBLE EULER LIMIT AND HILBERT'S EXPANSION 

Whenever a gas evolves in a fluid dynamic regime (at the length scale L), the average time between successive 
collisions involving a typical gas molecule is much smaller than the time necessary for an acoustic wave to travel 
a distance L. In other words, fluid dynamic regimes are characterized by the condition T c <C T a , or equivalently by he 
condition Kn <C 1 . 

In [22], Hilbert studied the Boltzmann equation (5) in the asymptotic regime defined by Kn = e <C 1 and Sh = 1. 
His idea was to seek the solution F £ of 

d t F e + v-V x F e = l^(F e ) (9) 



as a formal power series in e with smooth coefficients — known as Hilbert's expansion: 

F £ (t,x, v) = £ "fn{t,x, v) , with /„ smooth in t,x,v, foreachn>0. 

n>0 

He found that the leading order term in that expansion is of the form 

fo(t,x,v) =^( p ,„,e)( M: )(v), 
where (p,u, 9) is a solution of the compressible Euler system 

d t p+diw x (pu) =0, 

p(d,u + u-V x u)+V x (pd)=0, 

d t d + u-V x 6 + \6dxv x u = . 



(10) 



(11) 



Caflisch [10] succeeded in turning Hilbert's formal result into a rigorous statement bearing on solutions of the 
Boltzmann equation, by using a truncated variant of the Hilbert expansion above. Specifically, given a smooth solution 
(p , u, 9) of the compressible Euler system on some finite time interval [0, T), he constructs a family of solutions of the 
Boltzmann equation that converges to ^#( p , H] e) uniformly in t e [0, T) as e — > 0. Before Caflisch's result, Nishida [31] 
had proposed another proof of the compressible Euler limit of the Boltzmann equation under more stringent regularity 
assumptions, viz. analyticity, using some abstract variant of the Cauchy-Kowalewski theorem. 

One striking advantage of the Hilbert expansion is its versatility, abundantly illustrated by the great diversity 
of physically meaningful applications to be found in the work of Sone [36, 37]. However, there are some serious 
difficulties with the Hilbert expansion, some of which can be treated with adequate mathematical techniques. First, the 
radius of convergence of the Hilbert power series is in general, so that essentially all mathematical arguments based 
on Hilbert's expansion use a truncated variant thereof. In general, truncated Hilbert expansions are not everywhere 
nonnegative, and are not exact solutions of the Boltzmann equation. One obtains exact solutions of the Boltzmann 
equation by adding to the truncated Hilbert expansion some appropriate remainder term, satisfying a variant of the 



Boltzmann equation that becomes weakly nonlinear for small enough e (see for instance [10, 11, 2].) The truncated 
HUbert expansion with the remainder term so constructed is a rigorous, pointwise asymptotic expansion (meaning that 
£~ n \F e - (/o + e/i + ... +£ n fn)\ -> pointwise in (t,x,v)) of the solution F e of (9) as e -> 0. Another difficulty in 
working with Hilbert's expansion, even truncated at some finite order, is that /„ = 0(|V?„/o|) for each n > 0. Since 
generic solutions of the compressible Euler system lose regularity in finite time [33], truncated HUbert expansions 
make sense on finite time intervals only. For instance, if a solution (p , u , 9 ) of the compressible Euler system involves 
a shock wave, only the 0-th order term in the associated Hilbert expansion, i.e. fo(t,x,v) = •^ip, u ,8)(t,x)( v ) is weu 
defined. In general, if the geometric structure and the position of the singularities in the solution of the hydrodynamic 
equations are known precisely, one can bypass this difficulty by adding to the truncated Hilbert expansion appropriate 
boundary layer terms. If the structure of these singularities is unknown, or one does not even know whether the 
hydrodynamic solution is smooth, one cannot use the Hilbert expansion. 



GLOBAL EXISTENCE THEORY FOR THE BOLTZMANN EQUATION 

To avoid the various shortcomings of the Hilbert expansion method, one needs a theory of global solutions for the 
Boltzmann equation based on the only estimates that are uniform in Kn as Kn — > 0. These estimates are those deduced 
from the conservation laws (2) and Boltzmann's H Theorem (3), or from their differential formulations (7)-(8). 

Henceforth, we are concerned with solutions of the Boltzmann equation for a gas filling the Euclidian space R 3 and 
at equilibrium at infinity. By Galilean invariance and with a convenient choice of units, we can assume without loss 
of generality, that this equilibrium state at infinity is the Maxwellian ^#(1,0,1)) denoted by M in the sequel. In other 
words, we seek the solution of 

Shd,F + vV x F = J-^(F), (x,v) eR 3 xR 3 , t >0, 
Kn 

F{t,x,v)-*M as|x|^+°°, ( 12 ) 

F I — F in 
r \t=0 ~~ r ■ 

A convenient quantity measuring the distance between two distribution functions in the context of the Boltzmann 
equation is the relative entropy: for F = F(x,v) > and G = G(x, v) > a.e. in (x, v) g R 3 x R 3 , 

H(F\G):= [[ {Fki(F/G)-F + G){x,v)dxdv. (13) 

J JR 3 xR 3 

Notice that a\n(a/b) —a+b>0 for each a > and b > 0, with equality if and only if a = b. Hence the integrand is a 
nonnegative measurable function and H(F\G) —Oif and only if F = G a.e. on R 3 x R 3 . 

Since InM = — |ln(2jr) - 3M 2 , a formal argument based on the local conservation laws (7) and the differential 
inequality (8) shows that any classical solution F of (12) with appropriate decay as |v| — > +°° satisfies 

Shd, [ (F\n(F/M)-F+M)dv + dw x [ v(Fln(F/M) -F+M)dv < 0. 
Jr 3 Jr 3 

Integrating in x both sides of this inequality and assuming that F — >• M fast enough as |jc| — > +°°, we conclude that 

snpH(F(t,-,-)\M) <H{F in \M), and f ff -^(F)lnFdsdxdv <ShH(F in \M) . (14) 

t>Q JQ JJR 3 xR 3 

Observe that the collision integral ^(F) acts as a nonlocal integral operator analogous to a convolution in the v 
variable and as a pointwise product in the x variable. The fact that ^(F) is quadratic in F while H(F\M) is "essentially 
homogeneous of degree 1 as F 3> 1" suggests that c & (F) may not be defined for all nonnegative measurable functions 
F satisfying the entropy bound (14) above. Yet, for each measurable F > on R 3 x R 3 , one has 

^^[dxdv<C [ [ (-tf(F)lnF + (I + \v\ 2 )F)dxdv 



\x\ + \v\<r y/T+F J J\x\<r 



so that "io{F)/^\ + F e Lj 0C (R + x R 3 x R 3 ), i.e. is locally integrable in (f,x, v). This suggests dividing both sides of 
the Boltzmann equation by VI +F, thereby leading to the notion of renormalized solution. 

Definition. (DiPerna-Lions [12]) A renormalized solution relative to M of the Boltzmann equation is a nonnegative 
function F e C(R + ,L/ oe (R 3 x R 3 )) satisfying H(F(t)\M) < +°° for each t > and 

M(Sh d, + v • V X )F(F/M) = ^{F)T' (F/M) 

in the sense of distributions on R* + xR 3 x R 3 , for each T e C 1 (R+) satisfying T'(Z) < C/y/l+Z. 

With this notion of solution, one can prove the global existence and weak stability of solutions of the Cauchy 
problem for the Boltzmann equation, with initial data that are not necessarily small perturbations of either the vacuum 
state or of a Maxwellian equilibrium. 

Theorem. (DiPerna-Lions-Masmoudi [12, 27, 29]) For each measurable F m > a.e. on R 3 x R 3 satisfying the 
condition H(F'"\M) < +°°, there exists a renormalized solution of the Boltzmann equation (12) with initial data 
F m . This solution F satisfies 

Sh<9 ( / Fdv + div x [ vFdv = 0, 

, r JR3 (15) 

Shd, vFdv + div x v ® vF dv + div x m = , 
Jr 3 JR 3 

where m — m T > is a matrix-valued Radon measure, and the entropy inequality 

ShH(F(t,-,-)\M)+Sh f Tr(m(t)) - [' [ [ tf(F)\aF(s,x, v)dsdxdv < ShH(F in \M) , t>0. (16) 
Jr 3 JoJJr 3 xr 3 

A classical solution of the Boltzmann equation with appropriate decay as |v| — > +°° would satisfy all these properties 
with m~0; besides the entropy inequality is a weakened variant of Boltzmann's H Theorem — which would imply 
that this inequality is in fact an equality. 

The main advantage of the notion of renormalized solutions is that a) such solutions always exist for each initial 
data with finite relative entropy with respect to M, and b) such solutions are weakly stable, in the sense that if a 
sequence (F n ) n >o of renormalized solutions of the Boltzmann equation converges to F in the sense of distributions 
and satisfies H(F„\ t=Q \M) < C for all n > 0, where C is some positive constant, then F is also a renormalized solution 
of the Boltzmann equation, satisfying (15) and (16). Unfortunately, there is no uniqueness theorem for this notion of 
solution, so that a renormalized solution of the Boltzmann equation is not completely determined by its initial data. 
But if the Cauchy problem for the Boltzmann equation has a classical solution F, each renormalized solution of the 
Boltzmann equation with the same initial data as F coincides with F a.e. in (t,x,v) (see [26].) 



FLUID DYNAMIC LIMITS OF THE BOLTZMANN EQUATION 

As explained above, all fluid dynamic limits of the Boltzmann equation are characterized by the scaling condition 
Kn <C 1 : hence we set Kn = e throughout the present section. 

Besides, all the fluid dynamic limits considered in this paper correspond with weakly nonlinear regimes at the 
kinetic level — which does not imply that the nonlinearities are weak at the macroscopic level. Such regimes have 
been systematically explored by Sone at the formal level (see [36] and the references therein), by using the Hilbert 
expansion method. In other words, the distribution functions F considered are small perturbations of the Maxwellian 
state M at infinity. Henceforth, we denote by 8 £ <C 1 the order of magnitude of the difference F — M. A typical example 
of such a distribution function is F(t,x,v) — ^n,8 eU (t,x),i)( v )> since 

In this example, the distribution function F defines a velocity field u F and a temperature field 9f by the formulas 




Introducing the speed of sound Cf := ^/50f/3, we see that the Mach number Ma := Uf/cf = 8 e u, so that the scaling 
parameter 8 £ can be thought of as the (order of magnitude of the) Mach number. 



The acoustic limit 



The acoustic limit is the linearized variant of the compressible Euler limit considered by Hilbert himself. 

Theorem. (Golse, Jiang, Levermore, Masmoudi [14, 23]) Let Kn = £, Ma = 8 £ = 0(^/s) and Sh = 1. For each 
p' n ,u m ,9"' e L 2 (R 3 ), let F e be a family of renormalized solutions of the Boltzmann equation (12) with initial data 

K" = ^{i+s e p in ,s e u in ,\+s e e in ) ■ 

Then, in the limit as e — > 0, 

1 / ( J F £ (f,x,v)-M(v))(l,v,i|v| 2 -l)dv^(p, M ,0)(f,x) 
in Lj oc (R* + x R 3 ), where (p,u, 9) is the solution of the acoustic system 



d,p +div x M = 0, 


P\t=0 


= p in 


d,u + V x (p + 9)=0, 


U \ t =0 


= u in , 


d,d + \d3\ x u = Q, 




= 9 in 



While the result in [14] holds for the most general class of molecular interactions satisfying some angular cutoff 
assumption in the sense of Grad [20] (in fact, a much weaker version of Grad's assumption [23, 25]), an earlier 
contribution of the same authors with Bardos [5] introduced a key new idea in the derivation of hydrodynamic limits 
of renormalized solutions of the Boltzmann equation and treated the case of bounded collision kernels (e.g. cutoff 
Maxwell molecules). 



The incompressible Euler limit 

It is a well-known fact that, in the low Mach number limit, the flow of an inviscid fluid can be approximately 
decomposed into its acoustic and vortical modes, whose interaction vanishes with the Mach number. The result below 
explores the counterpart for vortical modes of the acoustic limit of the Boltzmann equation. Because of the low Mach 
number scaling, vortical modes evolve on a longer time scale than acoustic modes, consistently with the fact that the 
conditions V x (p + 9) — and div^M = characterize the equilibrium points of the acoustic system. 

Theorem. (Saint-Raymond [32]) Let Kn = e, and Sh = Ma = 8 e = e a with < a < 1. Let u in e H 3 (R 3 ) 1 satisfy 
diva" 1 = 0, and let u e C([0, T];H 3 (R 3 )) be a solution of the incompressible Euler equations 

J d t u + u ■ V x u + V x p = 0, div x M = 0, 

Let F e be a family of renormalized solutions of the Boltzmann equation (12) with initial data 

F e" = {l, Sett", l) ■ 

Then, in the limit as e — > 0, one has 

1 ( vF e (t,x,v)dv^u(t,x) inL~([0,r];LL(R 3 )). 

£ JR 3 

The proof of this result is based on the relative entropy method, described in the next section. Actually, there had 
been precursors of this theorem due to the author [9] and to Lions-Masmoudi [29], where the relative entropy method 
was introduced for this type of problem. Unfortunately, the statements in [9, 29] rested on extra assumptions on the 
family of solutions of the Boltzmann equation that remain unverified. These assumptions were removed by some clever 
argument in [32], which therefore contains the first complete proof of the theorem above. 



1 The notation H m (R" ) designates the Sobolev space of square integrable functions on R" whose partial derivatives of order < m in the sense of 
distributions are square integrable functions on R". A vector field is said to belong to H m (R") if all its components belong to H m (R"). 



The Stokes limit 



We continue our exploration of vortical modes with the Stokes limit of the Boltzmann equation. The scaling is 
weakly nonlinear at the macroscopic level of description, and the time scale is chosen so as to keep track of entropy 
production in the fluid dynamic limit. 

Theorem. (Golse, Levermore, Masmoudi [14, 25]) Let Kn = Sh = e, and Ma = 8 £ = o(e). For each (u in ,d in ) e 

L 2 x L°°(R 3 ) such that div x ii m = and each e e (0, 1 1 1/ 1 " 1 1 z.°° ), let F e be a family of renormalized solutions of the 
Boltzmann equation (12) with initial data 

Fg" = ^(i-8 e ein,5 e u in ,i+S e e in ) ■ 

Then, in the limit as e — > 0, one has 

l^ 3 ( J F £ (f,x,v)-M(v))(v,I|v| 2 -l)dv^(«,0)(f,x) inL; ot .(R + xR 3 xR 3 ), 

where (u, 0) is a solution of the Stokes-Fourier system 

f d t u + V x p = vA x u , diw x u = , u L_q — w m , 

The viscosity and heat conductivity in this theorem are given by the formulas (equivalent to the usual ones in [36]): 

v=i^*(v®v-i|v| 2 /), JC=§^*(i(|v| 2 -5)v), (17) 
where Si* denotes the Legendre dual of the Dirichlet form 3> of the collision operator linearized about M, i.e. 

0(<&):=i /// |4> + 4> H ,-4>'-4> / H ,| 2 |(v-v H ,)-co|MM < ,dvdv <I da). 

JiJR 3 xR 3 xS 2 

The fluid dynamic model obtained in the statement above is the Stokes-Fourier system; notice that the motion and 
temperature equations are decoupled in the absence of an external force field deriving from a potential. Previously 
Lions and Masmoudi [29] arrived at the particular case of the statement above corresponding to an initial data for 
which d m = 0, leading to the motion equation only, i.e. the evolution Stokes equation. For want of a better control 
of the high speed tails of the distribution function, their argument cannot be generalized to obtain the Stokes-Fourier 
system presented above. The proof in [14] uses a different idea originating from [5]. 



The incompressible Navier- Stokes limit 

Finally, we remove the weakly nonlinear scaling assumption at the macroscopic level of description, while keeping 
entropy production effects at leading order, and obtain the incompressible Navier-Stokes equations as a fluid dynamic 
limit of the Boltzmann equation. 

Theorem. (Golse, Saint-Raymond [18, 19]) Let Kn = Sh = Ma = 8 £ = £. For each (u in , 9 in ) e L 2 x L°°(R 3 ) such that 
div x u m = 0, let F £ be a family of renormalized solutions of the Boltzmann equation (12) with initial data 

F' £ n = . ^( i -eei»,£u' n ,\+£8' n ) ' 
for each e g (0, There exists at least one subsequence e„ — > such that 

l~ I (F e „(f,x,v)-M(v))(v,||v| 2 -l)dv^( M ,0)(f,x) inweak-4 e (R+xR 3 xR 3 ), 

W7 ^R 

where (u, 9) is a "Leray solution" of the Navier-Stokes-Fourier system with viscosity v and heat conductivity K given 
by formula (17): 

{d t u + dbt x {u®u) + V x p — vA x u , div x M = 0, u\ (=0 = u m , 

d t e+div x (ue) = lKA x e, e| r=0 = 0''". 



TABLE 1. Fluid dynamic limits of the Boltzmann equation, depending on 
the dimensionless parameters Kn, Ma and Sh. 



Boltzmann equation Kn = e <g 1 


| Ma i Sh 


Fluid dynamic limit 


| <5 e «:^/e | 1 


] Acoustic system 


) S £ < e | e 


! Stokes-Fourier system 


| 8 £ = e a , < a < 1 | 8 £ 


j Incompressible Euler equations 


j £ | £ 


| Incompressible Navier-Stokes equations 



Let us briefly recall the notion of Leray solution of the Navier-Stokes-Fourier system. In [24] (arguably one of the 
most important papers in the modern theory of partial differential equations), Leray defined a convenient notion of 
weak solution of the Navier-Stokes equations, and proved that, in space dimension 3, any initial velocity field with 
finite kinetic energy launches at least one such solution defined for all times. Leray solutions are not known to be 
uniquely defined by their initial data; however, if an initial data launches a smooth solution, all Leray solutions with 
the same initial data must coincide with that smooth solution. At the time of this writing, it is yet unknown (and a 
major open problem in the analysis of partial differential equations) whether Leray solutions launched by any smooth 
initial data remain smooth for all times. Thus, we do not know whether different subsequences £„ — > in the theorem 
above lead to the same Leray solution (u, 9) of the Navier-Stokes-Fourier system in general. 

A Leray solution of the Navier-Stokes-Fourier system above is a pair (u, 9) consisting of a velocity field u and 
a temperature field 9, both continuous on R+ with values in L 2 (R 3 ) equipped with its weak topology, that solves 
the Navier-Stokes-Fourier system in the sense of distributions, satisfies the initial condition, and verifies the Leray 
inequality: 

\\ (| M | 2 + ||0| 2 )(f,x)dx+ f I (v\V x u\ 2 + K\V x 6\ 2 )(s,x)dxds<± [ {\u in \ 2 + l\6 in \ 2 ){t,x)dx. (18) 
Jr 3 Jo Jr 3 Jr 3 

The Leray inequality is an equality for classical solutions of the Navier-Stokes equations, exactly as the DiPerna-Lions 
entropy inequality (16) is an equality for classical solutions of the Boltzmann equation. This indicates that the Leray 
existence theory for the Navier-Stokes equations and the DiPerna-Lions existence theory for the Boltzmann equation 
are parallel theories. The theorem above explains how these theories are related in the hydrodynamic limit. 

Partial results on this theorem have been obtained by Lions-Masmoudi [28]. While the reference [18] treated the 
case of bounded collision kernels, the theorem above was later extended to all hard cutoff potentials in the sense of 
Grad — which includes the case of hard spheres considered in this paper — in [19]. The arguments in [18, 19] have 
been recently refined by Levermore and Masmoudi [25] to encompass both soft as well as hard potentials, under a 
cutoff assumption more general than that proposed by Grad in [20]. 

While these results bear on the most general case of renormalized solutions without restrictions on the size of initial 
data in space dimension 3, the Navier-Stokes limit of the Boltzmann equation had previously been obtained in the case 
of global smooth solutions for small initial data by Bardos and Ukai [7]. The Navier-Stokes limit of the Boltzmann 
equation had also been established on finite time intervals by adapting the Caflisch method based on Hilbert truncated 
expansions, by DeMasi, Esposito and Lebowitz [11]. 

The fluid dynamic limits discussed in this section can therefore be summarized as in table 1 . Notice that these limits 
have been established for molecular interactions more general than hard sphere collisions; see the references listed in 
the statements of the various theorems above for the conditions on the collision kernel b(v — v*, <X>). All these results 
assume some angular cutoff on the collision kernel as proposed by Grad [20] — or slightly more general, as in [25]. 

More importantly, some of the conditions bearing on the parameters Kn, Ma and Sh may be not optimal. Formal 
arguments suggest that the acoustic limit should hold whenever 8 £ <C 1 instead of 8 £ <C \fe, while the incompressible 
Euler limit should hold under the weaker condition 8 £ ^ £ instead of 8 e = E a with < a < 1. 

Let us conclude this section with an important remark on the physical meaning of the "incompressible" fluid 
dynamic limits of the Boltzmann equation. What is proved in the last three theorems is that, to leading order, the 
velocity field u satisfies the same equations as the velocity field of an incompressible fluid with constant density. This 
does not mean that the gas is incompressible in that regime. Also, in the case of an incompressible fluid with the same 
heat capacity and heat conductivity as the gas, the diffusion term in the equation for the temperature field would be 



multiplied by 5/3. This difference comes from the work of the pressure: see the discussion in footnotes 6 on p. 93 in 
[36] and 43 on p. 107 of [37], together with section 3.7.2 in [37]. 

Likewise, the inequality (18) was written by Leray in [24] with = 0. For an incompressible fluid with constant 
density p , the quantity j J p | u (t , x) | 2 dx is the kinetic energy of the fluid at time t , and the Leray inequality is interpreted 
as a statement on the dissipation of energy in the fluid. The meaning of (18) with 9 ^ is obviously different, since 
the quantity | / (|m(?,jc)| 2 + j9(t,x) 2 )dx is not the total energy of the gas at time t. 

MATHEMATICAL TOOLS FOR THE HYDRODYNAMIC LIMIT 
The linearized collision integral 

In all the fluid dynamic limits considered in the previous section, the solution F £ of the Boltzmann equation (12) 
is a small perturbation of the uniform Maxwellian equilibrium state M. Therefore, the linearization about M of the 
Boltzmann collision integral plays an important role in these limits. Thus, we consider this linearized collision integral 
intertwined with the multiplication by M, and set Jz4f = — M~ x D c <o (M) ■ (M0 ), or equivalently 

&M<l>(v) := [ [ (<j)(v) + <p(v,)-^(v')-<j)(v',))\{v-v,)-(o\M(v,)dv,d(0. (19) 

J JR 3 xS 2 

Hilbert [22] proved that «2m is an unbounded, Fredholm, self-adjoint nonnegative operator on L 2 (R 3 ;Mdv) 2 , with 
domain L 2 (R 3 ;(1 + |v| 2 )Mdv) and nullspace KerJz?M = Span{l, vi, V2, V3, |v| 2 }. Hilbert's argument, written for the 
hard sphere case, was later extended by Grad [20], who defined some appropriate class of collision kernels b(v — v* , a>) 
for which the linearized collision integral satisfies the Fredholm alternative. Grad's idea was that grazing collisions 
between neutral gas molecules are rare events that can be somehow neglected, at variance with the case of plasmas or 
ionized gases. Henceforth, we denote 

(0)=/ <j)(v)M(v)dv for each e L 1 (R 3 ;Mdv) . 

With this notation, the Fredholm alternative for the integral equation = S with unknown / and source term 
can be stated as follows: 

a) either (S) = (Svi) = (Sv 2 ) = (5v 3 ) = (S|v| 2 ) = 0, in which case the integral equation has a unique solution / 
satisfying 

/GL 2 (R 3 ;(l + |v| 2 )Mdv) and (/) = (/v x ) = (fv 2 ) = (/v 3 ) = (/|v| 2 > = 0, 

henceforth denoted / = Jz^'S, or 

b) there exists <E KerJz^ such that (50) 7^ 0, in which case the integral equation Jz?m/ = S does not have any 
solution in L 2 (R 3 ; (1 + |v| 2 )Mdv). 



The moment method for the Navier-Stokes limit: formal argument 

Define g e by the formula F e = M(l + 8 £ g £ ). If F £ satisfies (12) with Kn = Sh = 8 £ = £, the relative fluctuation of 
distribution function g £ satisfies 

ed tge + v ■ V xg£ + j&Mge - £M(ge,ge) , (20) 

where £} M is the symmetric bilinear operator defined by =2^(0,0) =M~ 1<r #'(M0). Multiplying each side of (20) by e 
and letting e — > shows that, if g £n — > g for some subsequence e„ —¥ 0, the limiting fluctuation g is a "local Maxwellian 
state", i.e. is of the form 

g(t,x,v) =p(f,x)+ M (f,x)-v+0(f,x)|(|v| 2 -3). (21) 



2 The notation If(R N ;fdv) (where p > 1 and / > is a measurable function defined a.e. on Q.) designates the set of measurable functions <f> 
defined a.e. on that satisfy 

/ |</>(v)F7(v)dv<+°°. 



Multiplying each side of (20) by M and vM, and integrating in v e R 3 shows that 

ed,{g e )+drv x {vg e ) = 0, ed,{vg e ) +div x {v®vg e ) = 

in view of (2). Passing to the limit as £„ — > 0, and taking into account the local Maxwellian form (21) of g leads to 

div x u = div x (vg)=0, V x (p + 9)=div x (v®vg}=0. (22) 

The first equality is the solenoidal condition for the velocity field u, while the second implies that p + = 0, assuming 
that p , 9 -> as |jc| -> +°°. 

Next we multiply each side of (20) by jvM and integrate in v e R 3 to obtain 

1 1 /Ivl 2 \ 

d t {vg e )+6iv x -{Ag e ) = -V x -(^j-g e ), 

where A(v) := v- ||v| 2 /. One has (A«) = (A«vi) = (A w v 2 ) = (A H v 3 ) = (A«|v| 2 ) = for each fc,/ = 1,2,3, so 
that A w := Jf M 'A w e L 2 (R 3 ; (1 + |v| 2 )Mdv) is well-defined. Since «5f M is self-adjoint on L 2 (R 3 ;Mdv), one has 



\{Ag e ) = \{{J? M A)g e ) = {A^ M8e j = (A(<2 M (g e ,g e )) - (A(ed t g e + vV xge )) ^ 

^(A(J2 M (g,g))-(AvV x g). 

By (21) and the solenoidal condition in (22), the second term takes the form 

(Av • V x g) = (A ® v ® v) : V xM = v(V xM + (V lM ) r ) . (24) 

Indeed (AyAj/) = v(8n c 8ji + 5 ( ;5^ — §5y 5^;), which can be recast as v = ^, (A : A) since A(/?v) = 7M(v)^ r for each 
v e R 3 and each R E Oj, (R). This formula for v is equivalent to the first relation in (17). 

As for the first term, since g € KerJ£ M according to (21), one has £} M (g,g) = \^£ M {g 2 ) (see [3], fla. (60) on p. 338.) 
Hence 

(A(£ M (g,g)) = \(A^ M (g 2 )) = \{{X u A)g 2 ) = \{Ag 2 ) = \(A®v®v) : u® u = u®u- ||«| 2 /, (25) 

in view of the elementary identity (AyA w ) = (8 ik 8ji + 8 a 8 jk - f 5 y 5 fc; ). 

Let E, = E, (x) e C~(R 3 ) be a divergence-free test vector field. Substituting (24) and (25) in (23) shows that 

0=d t Jt J - (vg e )6x-J Vi; : ^(Ag^dx 

-td, j %-udx- J VE, : {u®u-\\u\ 2 I)Ax + v J VE, : (V x u + (V x ii) T )dx. 
Since E, is divergence-free 

J VE, : \\u\ 2 I6x= \ J | M | 2 div x <^dx = 0, 

while 

J VE, : {V x ufdx=-J V(div<i;)- M dx = 0. 

Therefore 

d, J E, u6x- J VE, : u®udx + v J VE, : V x udx = (26) 

for each divergence-free test vector field E, = 4(x) e C~(R 3 ). Now, if T e £^'(R 3 ) is a vector- valued distribution 
satisfying (T,%) =0 for each divergence-free test vector field E, = E,(x) e C"(R 3 ), there exists a scalar distribution 
71 € f^'(R 3 ) such that T = Vtt. In other words, (26) means precisely that m is a weak solution of the motion equation 
in the Navier-Stokes system. 



Compactness tools 



An important ingredient in the proof of all fluid dynamic limits of the Boltzmann equation considered above is the 
fact that the relative fluctuation of distribution function g £ := (F £ —M)/8 £ M converges in some sense, possibly after 
extracting some subsequence e„ — > 0. The key argument is the following inequality resulting from (16): 

/ (h(8 ege )){t)dx = H(F e (t,;-)\M)<H{F™\M) = 0{8 2 e ) foreachf>0, 

for the initial data considered in the four theorems stated in the previous section, where h(z) := (1 +z)ln(l + z) — z. 
Since h(z) ~ z 2 /2 as z — > 0, this control is as good as a bound in L°°(R + ;L 2 (R 3 x R 3 ;Mdxdv)) for the values of 
g e not exceeding 0(1/ 8 £ ). Thus (1 + |v| 2 )g £ is relatively compact in weak-L 1 ([0, T] x [— R,R] 3 x R 3 ;Mdf dxdv) for 
each R,T > 0, and all its limit points as e — > belong to L°°(R + ;L 2 (R 3 x R 3 ;Mdxdv)). In the case of the acoustic 
or Stokes-Fourier limit, the uniqueness of the solution of the limiting fluid equations implies that the whole family 
(1 + |v| 2 )g e converges weakly. 

Since the leading order term in (20) is ^J^Mgs and _2m is a linear operator, the weak compactness of the family 
(1 + |v| 2 )g £ is enough to conclude that any limit point g of that family as£->0 must satisfy J^m8 — 0, and therefore 
is an infinitesimal Maxwellian, i.e. is of the form (21). 

In addition, for the Navier-Stokes-Fourier limit, the compactness of the family g £ in Lj oc for the strong topology 
(implying the a.e. pointwise convergence of a subsequence) is needed to pass to the limit in nonlinear terms. We 
use repeatedly some compactness results for moments of the distribution function in the velocity variable v based on 
bounds on the streaming operator — see [16, 15]. These compactness results are referred to as compactness by velocity 
averaging. A typical example of a velocity averaging theorem used in the Navier-Stokes limit of the Boltzmann 
equation is as follows. We state it in the steady case for the sake of simplicity. 

Theorem. (Golse, Saint-Raymond [17]) Let /„ = /„(x,v) be a bounded sequence in i ! (R^ x R^) such that the 
sequence v ■ V x f„ is bounded in L 1 (R^ x R^), while /„ itself is bounded in L 1 (R^ ;L P (R^)) for some p> I. Then 

a) /„ is weakly relatively compact in L l loc (R^ x R^); and 

b) for each <j) e C C (R N ), the sequence of velocity averages 

/ fn(x,v)<j)(v)dv 

JR N 

is strongly relatively compact in Lj oc (R N ). 



The conservation laws 

The formal argument presented above in the case of the Navier-Stokes limit shows the importance of the local con- 
servation laws of mass, momentum and energy in the derivation of fluid dynamic models from the Boltzmann equation. 
Unfortunately, renormalized solutions of the Boltzmann equation are not known to satisfy the local conservation laws 
of momentum and energy in (7). They satisfy instead the approximate conservation laws 




(27) 

Therefore, one must show that the conservation defects 

in L ; 1 oc (R + x R 3 ) as e — > 0, and identify the limits as e — > of the terms 



This raises an important question regarding the tail of the distribution functions F £ as |v| — > +°°. That the family 
(1 + |v| 2 )g £ is relatively compact in weak-L 1 ([0, T] x [— R,R} 3 x R 3 ;Mdf dxdv) for each R,T > is in general not 
enough — for instance, in the acoustic limit, Sh = 1 and one needs to identify the limit of the energy flux 

which is a 3rd order moment in v of ^-(r(l + 8 £ g £ ) — T(l)) <~ T'(l)g £ . Controlling the high speed tail of (fluctuations 
of) the distribution function is an essential step in the derivation of fluid dynamic limits of the Boltzmann equation, 
and involves rather technical estimates based on the entropy and entropy production estimates (16) together with the 
dispersion effects of the streaming operator Sh<9 r + v • V x (see [4, 18, 19]). 



The relative entropy method 

In inviscid hydrodynamic limits, i.e. the compressible or incompressible Euler limits, entropy production does not 
balance streaming. Therefore the velocity averaging method fails for such limits. The idea is to use the regularity 
of the solution of the target equations, together with relaxation towards local equilibrium in order to obtain some 
compactness on fluctuations of the distribution function. 

Pick for instance u, a smooth solution of the target equations — e.g. the incompressible Euler equations — and 
study the evolution of the quantity 

Z£ W = ^^(-^ 1*^(1 A«M,i)) 

where F £ is a renormalized solution of (12) with Sh = Ma = 8 £ = E a and < a < 1. This is the leading order 
of the relative entropy of the Boltzmann solution with respect to the local Maxwellian state defined by u, in the 
incompressible Euler scaling. At the formal level, it is found that 

^W = -i / V x ii: f {v-8 e u)® 2 F e dvdx+±- [ V x p- f (v- 8 e u)F e dvdx. 

at Og JR3 JR3 £ JR3 JR3 

The second term in the r.h.s. of the equality above vanishes with e since 

-j- f vF e (t,x, v)dv — > divergence-free vector field. 

8 £ JR3 

The key idea is to estimate the first term in the r.h.s. as follows 

ff \V x it : (v- 8 e uf 2 F e \ dvdx< CZ £ (t) +o(l) 

where C — 0(\\V x u\\l^). Then, one concludes with Gronwall's inequality. 

The relative entropy method stems from an idea of H.T Yau (for Ginzburg-Landau lattice models, see [39]); it was 
later adapted to the Boltzmann equation by the author [9] and Lions-Masmoudi [29]. It is especially designed to handle 
sequences of weak solutions of the Boltzmann equation converging to a classical solution of the fluid equation. 



CONCLUSIONS 

The DiPerna-Lions theory of renormalized solutions of the Boltzmann equation allows one to obtain derivations of 
fluid dynamic regimes from the kinetic theory of gases without unphysical assumptions on the size or regularity of the 
data. Following the program outlined in [4], these derivations are based on 

a) relative entropy and entropy production estimates, together with 

b) functional analytic methods in Lebesgue (L p ) spaces. 

At present, the program in [4] leaves aside the compressible Euler limit of the Boltzmann equation, or the asymptotic 
regime leading to the compressible Navier-Stokes equations. Little progress has been made on these issues since the 
work of Nishida [31] and Caflisch [10]. 



The problem of deriving fluid dynamic limits from the Boltzmann equation in the steady regime is also of consid- 
erable importance for practical applications. Formal results are of course well understood with the classical Hilbert 
or Chapman-Enskog expansion techniques — see the book of Sone [36]. Unfortunately, the theory of the steady 
Boltzmann equation with large data is not as mature as its counterpart for the Cauchy problem, in spite of interesting 
contributions by Arkeryd and Nouri [1], and there is no analogue of the DiPerna-Lions theory for the steady case yet. 

But even for evolution problems in regimes that are weakly nonlinear at the kinetic level, the relative entropy is not 
the solution to all difficulties. In several asymptotic regimes of the Boltzmann equation, the leading order and next to 
leading order fluctuations of the distribution function may interact to produce highly nontrivial macroscopic effects in 
the fluid dynamic limit. Examples of such asymptotic regimes are 

a) ghost effects, introduced by Sone, Aoki, Takata, Sugimoto and Bobylev in [34], reported in Sone's Harold Grad 
Lecture [35] and [36, 37], and systematically studied by Sone, Aoki and the Kyoto school, 

b) Navier-Stokes limits recovering viscous heating terms, due to Bobylev [8] and Bardos-Levermore-Ukai-Yang [6] 
— see also the discussion in [38], and 

c) hydrodynamic limits for thin layers of fluid — see [13]. 

Acknowledgements. The author thanks Profs. Aoki, Levermore and Sone for their generous scientific advice during 
the preparation of this paper. 
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